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Inversion of Radio Wave Absorption Data To Establish 
Ionospheric Properties. I. Nondeviative Absorption 

Albert D. Wheelon 

(May 15, 1963) 

Nondeviative absorption of high frequency radio waves which penetrate the lower 
ionospheric layers is analyzed in terms of the usual expression : 
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Classical approximate predictions of the absorption function A(co) for several ionospheric 
electron density N and collision frequency v models derived from this expression are first 
reviewed. It is then shown that by solving t he above integral equation, one can establish an 
explicit relationship between N and v at any height from a measured knowledge of the 
absorption function A(a>) as a function of carrier frequency o>. Thus, a knowledge of the 
electron density profile can be used to reconstruct the collision frequency profile and con- 
versely. It is shown that the usual assumption A(u>) = Aqco' 2 corresponds to an impossible 
ionospheric model. Several examples are given to illustrate the analytical inversion tech- 
nique and to indicate its sensitivity to low frequency portions of the measured absorption 
function. 



1. Introduction 

The absorption of radio waves reflected from the 
ionosphere has been important for several decades 
in connection with estimating the required power 
margins of shortwave communication links. More 
recently, it has been recognized that the variation of 
radio wave absorption with carrier frequency and 
other independent propagation variables represents 
a valuable new tool for studying the electron density 
and/or collision frequency profile of the ionosphere 
[Little, 1961]. The general absorption problem 
divides itself naturally into two topics: deviative and 
nondeviative. The process is called deviative if the 
wave is reflected by the layer which is primarily 
responsible for the absorption. Deviative absorp- 
tion is the more difficult of the two types, since one 
must then recognize both the absorption and refrac- 
tive effects simultaneously at each point along the 
ray trajectory. Nondeviative absorption refers to 
wave absorption by layers which lie well below the 
reflecting layers ; in this case the ray paths are neither 
bent nor "stretched" in the region of principal 
absorption. Of course, a radio star or satellite signal 
which passes completely through the ionosphere is 
also beset by nondeviative absorption. Figure 1 
illustrates these three cases. 

Tt is commonly held that the D region of the iono- 
sphere between 70 and (SO km is primarily responsible 
for shortwave absorption. The absorption per unit 
length is roughly proportional to the product of 
electron density AT" and collision frequency v, and the 
D region is that in which the suddenly rising electron 
density profile and exponentially declining collision 
frequency have the greatest product. For vertical 



incidence, the deviative-to-nondeviative frequency 
"watershed" occurs at several megacycles, corre- 
sponding to the critical frequency of the E region 
l>\ day. 

This first paper is concerned with nondeviative 
absorption, and is therefore relevant to both the 
absorption of radio star /satellite signals propagating 
through the ionosphere and F layer reflection by day 
and night. We base our treatment here on the 
Appleton-Hartree expressions for the electromagnetic 
response of an ionospheric plasma which suffers 
collisions with the neutral particles in a way which 
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Figure 1. Typical ray paths for radiowave absorption for rays 
entering the absorbing D-region with an angle of incidence i: 
(a) deviative absorption in reflection region; (b) nondeviative 
absorption below reflecting F layers, and (c) nondeviative 
absorption of a radio star or satellite signal. 
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does not depend on the energy of the electrons them- 
selves. We shall return in a later paper to consider 
the effect of modifications of this analysis called out 
by the revision of the Appleton-Hartree formula 
which recognizes an electron energy dependence of 
the collision cross section with neutral atoms [Sen 
and Wyller, I960]. 

The starting point for this analysis is the magneto- 
ionic expression for field strength attenuation [Rat- 
cliffe, 1959]. Above several megacycles, it is probably 
safe to use the quasi-longitudinal approximation, 
which relates the received field strength E R (cc) to 
that which would be received if there were no 
absorption E (o)) by the following expression: 



E R (a>)=E (o>) exp- 
where 
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e=4.8X10 -10 esu, electronic charge. 
m = 9.11X10 -28 g, electronic mass. 
c = 3X10 10 cm/sec, speed of light. 
co L = Local gyrof requency of the earth's magnetic 
field times the cosine of the angle \p between 
the local ray direction and the magnetic 
field vector. 
v= Local collision frequency. 
N= Local electron density. 
id = Corresponding local refractive index. 

The integral in (1.1) is to be taken along the nominal 
ray path determined by SnelPs law and the refractive 
index profile n(y). For the nondeviative case under 
consideration, the refractive index may be set equal 
to unity, and the integral along the ray path in the 
absorbing region replaced by that along a rectilinear 
path, 

ds = dy sec i 

where i is the angle of inclination of the ray with re- 
spect to the local vertical (y) as it traverses the ab- 
sorbing region (see fig. 1). This allows (1.1) to be 
rewritten as: 



E B (Q})=E (a>) exp- 
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The height range of integration here has been ex- 
tended from ground level to infinity, since the non- 
deviative absorption is confined to the lower reaches 
of the ionosphere and negligible error is introduced 
by extending the integral indefinitely above this. 
Strictly speaking, expression (1.2) corresponds to the 
radio star/satellite case shown in figure 1 and one 
should insert an additional factor of two in the ex- 
ponent to account for the double ionospheric traverse 
associated with F layer reflection. 

2. Classical Predictions of Absorption 

Classical predictive analysis of nondeviative ab- 
sorption approximates the integral term in (1.2) by 



assuming that the carrier frequency of the wave is 
considerably greater than the collision frequency 
over the height range for which the product Nv is 
significantly large. The same assumption implies 
that co^coc, and for F layer (reflective) propa- 
gation one obtains: 

E R (u)=E,{u) exp-(^) sec i — f °° dyNv. (2.1) 

Since the secant of the angle of incidence is roughly 
proportional to the transmission distance D for flat 
earth propagation via thin layers, one obtains the 
usual absorption expression [Best and Ratcliffe, 
1938]. 

E R =E exp-rjDX 2 . (2.2) 

This simple law has been confirmed approximately 
by numerous shortwave experiments [Farmer and 
Ratcliffe, 1935]. It is precisely this scaling law which 
encourages one to work as close as possible to the 
maximum usable frequency (MUF) on any given 
day, so as to reduce the effects of nondeviative 
absorption. 

A significant diurnal variation of the measured 
exponential coefficient 77 in (2.2) is properly attrib- 
uted to the corresponding diurnal variability of 
the electron density N in the defining integral (2.1). 
This variability was calculated theoretically by 
Appleton [1937] using the Chapman model for the 
electron density versus height profile [Chapman, 
1931]. 

N(y, X) =A^ exp-i QL + a Po H sec xe~^ (2.3) 

where X is the time-variable solar zenith angle, H 
the scale height of the exponential atmosphere 
assumed, A the photo-absorption coefficient, and p 
the sea level neutral particle density. Appleton 
evaluated the diurnal variability of the absorption 
factor by combining (2.3) with an exponential col- 
lision frequency profile 1 p(y) = v e~ v/H . Carrying 
out the integral in (2.1) Appleton [6] found 



E R c^E exp—rj'D cos 3/2 xX 2 . 



(2.4) 



1 Which is also implicit in the Chapman model of (2.3). 

2 Ionospheric radio propagation (June 1948), NBS Circ. 462, 37-38; 112-113. 



A good deal of effort has recently been devoted to 
comparing accumulating experimental data with 
expression (2.4). This effort has been moderately 
successful and many workers have been content to 
let the matter rest there. On the other hand, the 
extensive NBS empirical studies favor a linear varia- 
tion of the absorption factor 2 with cos X, rather 
than the cos 3 / 2 x term in (2.4). Beynon and Davies 
[1954] found that a frequency and zenith angle de- 
pendent coefficient of X 2 , A+BF(j) cos x, gives a 
better overall fit to the experimental data gathered 
on 2.0 Mc/s in England. Davies and Hagg [1955] 
have even challenged the X 2 scaling law of (2.4) after 
careful examination of absorption data on several 



700 



frequencies at high latitudes (Prince Rupert), with 
corresponding variation in the exponent in cos w x 
between m = 0.29 and 0.82. Chapman and Davies 
[1958] have also commented on the approximate 
constancy of daytime absorption. 

This variability and/or uncertainty probably 
represents a fundamental deficiency in the simplified 
prediction (2.4). This deficiency can be traced, at 
least in part, to the unwarranted simplification of the 
original expression (1.2). The assumption that "the 
collision frequency v is small compared with the 
carrier frequency co in those regions where Nv is 
large," is a qualitative argument whose precision 
and value has lost ground as more accurate models of 
the lower ionosphere have become available. 
Another way of looking at the matter of predominant 
absorption heights is to note that the factor 
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has its maximum value at the height corresponding 

to v=oo, and so long as iV has a respectable value at 
such heights (i.e., the D region for HF), the fre- 
quency variable denominator in (1.2) ought not to 
be approximated to the form of (2.1). In point of 
fact, it is just this denominator — in a delicate 
balance — which can account for the spectrum of 
frequency and zenith angle scaling law exponents 
observed experimentally. 

Jaeger [1947] published an important analysis of 
nondeviative absorption based on the complete ab- 
sorption integral in (1.2) which demonstrates this 
natural variability. His calculation is based on the 
Chapman electron density model (2.3) and an ex- 
ponential collision frequency profile, 

E B ~E exp-[„Z? cos 3 ' 2 X X 2 J (°^~j] (2-5) 
where 
J03)=4tt/3 3 "V^ {cos /3 [|- 0(0)] 

+sin/?Q-S(0)]| (2.6) 
m \ If* 7, COS t 
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The function J($) ranges from to 4.133 as varies 
from zero to infinity. The small argument behavior 
of this function, 

lim J(p)m2wy[ep /2 
indicates that the absorption exponent is independent 



of solar zenith angle (i.e., time) for small /3. This 
behavior thus provides a flexible range of frequency 
and zenith angle scaling law exponents, 

\ 1/2 D to \ 2 Z?cos 3/2 X 

as suggested by experiment. However, it does not 
yet seem that this broader theory has been put to the 
test of close comparison with experimental data. 

The real point of the matter would seem to be that 
theory has reached the point of diminishing returns 
in fitting theoretical electron density and collision 
frequency profiles to experimental absorption data, 
either through the complete magnetoionic expression 
(2.4) or its approximate form (2.1). The present 
troubles run sufficiently deep to encourage one to 
take another approach to the entire problem and 
provide the mot i vat ion of the present series of papers. 

In this first paper we shall regard absorption 
measurements as a means of establishing the electron 
density profile, rather than as a means of verifying 
preconceived notions of how this profile ought to 
look. In this sense, we shall regard (1.2) as an 
integral equation for finding A r (//) from a measured 
know ledge of E R /E Q as a function of carrier frequency. 
We shall show that this integral equation can he solved 
explicitly, thereby circumventing the need for fitting 
adjustable parameters in theoretical profile models. 

3. Inversion of the Absorption Integral 

Let us now turn to the direct problem of recovering 
the electron density or collision frequency profile 
from a measured variation of absorption with carrier 
frequency. Actually, we shall find a relation be- 
tween electron density and collision frequency at any 
height as the solution of the integral equation: 
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(3.2) 



is related to the basic measurements of field strength. 
Let us assume that the function A(u>) is measured 
by carefully calibrated experiments and show how 
one can solve the integral (3.1). 

The first step is to convert the integral in (3.1) 
over all height ranges (0<C^/<C°°) to an integral over 
the collision frequency v between corresponding 
limits. The value of v for an infinite height above 
the earth is zero, since the atmosphere constituents 
with which the electrons make collisions decay 
rapidly with altitude. Even in the case of F layer 
propagation (fig. lb), the collision frequency at the 
top of the ray path is probably 10 3 c/s, and therefore 
several orders of magnitude smaller than the mega- 
cycle values for v found in the absorbing (D) region. 
The ground level value of collision frequency is 
very much greater than any carrier frequency (i.e., 
10 23 ) and is therefore reasonably taken to be infinity. 
These limit considerations are crucial to what fol- 
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lows, but appear to be amply justified by the physics 
of the ionosphere as we know it today. With this 
change of independent variable, one is able to write, 



A(aj)== Jo rf V+'-- 1 -'- v 



-(uicoj 2 



(3.3) 



Where the derivative of v with respect to height is 
denoted by: 

, _dv 

and N is to be regarded as a function of collision 
frequency until the solution is effected. It seems 
that this trick requires that iV be a monotonic func- 
tion of v, although closer analysis may supply a 
more generous restriction. We next note the 
following (Laplace transform) identity: 



wicoz, 



v 2 +(co±co L ) 2 



■;; 



drje^ u±a ^ cos (rjv). (3.4) 



If this representation is substituted into (3.3) and 
the order of integration reversed, one finds 

/» CO /» CO 

(a)±u L )A(a))= dv 4 N{y) drje'^^^ COS (rjv) 
Jo v Jo 

drje-^Ae^^ I dv cos (tjv) ^N{v) 1- (3.5) 

This reversal indicates that the unknown function of 
7j in square brackets is the Laplace transform of the 
known function (o>±co L )A(u) on the left hand 
side of the equation. Inverting this relation, we find 
that 

e**»i<JC- 1 [(w±u L )A(u)]= f dv cos (vrj) ^N(v) t 

Jo v 

(3.6) 

where J0~ l [/] indicates the inverse Laplace transform 
of /. Equation (3.6) states that the Fourier cosine 
transform of the unknown function v/v'N(v) is equal 
to the known function on the left-hand side. Using 
the Fourier orthogonality result, 



I 



drj cos (art) cos (br))=ir5(a—b) (3.7) 



where 5(x) is the singular Dirac delta function, one 
finds that 



V T Jo 



dxcos (xv)e± x "L ( £- 1 [( 0) ±a) L )A(Go)] x , 



where the inverse Laplace transform is to be ex- 
pressed as a function of the dummy variable of 
integration x. Separating terms, we have the final 
solution 

N ^=- ( i- P<tecos (xv)e±*«L£-\(o>±o> L )A{o>)] x . 

(3.8) 



This equation provides an explicit relationship 
between the value of electron density N at an arbi- 
trary height y and the collision frequency v at the 
same height, and is therefore the desired solution. 
If the collision frequency profile is known as a func- 
tion of height, this expression gives an explicit elec- 
tron density profile. On the other hand, if it is the 
electron density profile which is known, one can 
reconstruct the collision frequency variation with 
height by integrating both sides of (3.8) from corre- 
sponding values of y and v 

/»CO S* v (y) j j /»CO 

dyN{y) = — dv-- \ dx 

Jy Jo J'TrJo 

cos {vx)e± x «L £-\(<*±<a L )A{u)] x . (3.9) 

This result completes the formal solution to the 
problem at hand. 

4. Illustrative Examples 

To show how the formal solution (3.8) can be used 
in a practical way to invert experimental measure- 
ments of the absorption function A(a>), it is instruc- 
tive to consider several examples, in which we shall 
take w L =0 for convenience. These examples will 
also indicate the extraordinary sensitivity of the 
electron density profile to the low frequency behavior 
of the absorption function. 

A. Example 1. Let us consider first the logical 
consequences of the simplified classical absorption 
result [cf. (2.2)]. 



A(«)=4 



O co 



(4.1) 



which is plotted in figure 2. Noting that the inverse 
Laplace transform of a? -1 is just unity, we have 

=A Q - — dx cos (vx) 1 

Jo 



' dy -k , 



or using (3.7), we have: 



N(y)-A- v ^Kv-o). 



(4.2) 



This result insists that the electron density be rigor- 
ously zero everywhere except at those (infinite) 
heights for which the collision frequency vanishes 
completely, where it must be infinite (see fig. 3a). 
Such a result is both ridiculous and inescapable, 
since it follows inevitably from the assumed fre- 
quency dependence (4.1) via the exact solution (3.8) 
of the absorption integral (3.1) describing the proc- 
ess. One can verify that (4.2) is the solution of 
(3.1) by substituting it therein and easily recovering 
the result (4.1). 

The fault lies with the implicit assumption that(4.1) 
describes the frequency variation of the absorption 
factor for all frequencies. It certainly does not 
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Figure 2. Plots of various absorption functions all of which 
have the same high frequency behavior. co~ 2 

apply to the low frequency end of the carrier spec- 
trum where this behavior must eventually change 
as one approaches the critical frequency of the 
absorbing layers, since below this plasma cutoff, 
one has no transmission to be absorbed. In point 
of fact, it is the high frequency asymptotic behavior 
of the absorption integral (3.1) which is expected to 
correspond to the simplified classical form (4.1). 
We therefore look for other models of the absorption 
factor A (od) which vary as co~ 2 for large values of the 
argument but which have more reasonable vari- 
ations for small w. 

B. Example 2. The following function removes 
the low frequency divergence of example (4.1), but 
in a particular way which controls the solution. 



A((a) = 



Aa 



~* 2 +tf 



(4.3) 



Taking this together with (3.4) gives 



-A i — -• f 



dx cos (xv) cos (xp) 




Figure 3. N versus v relationships corresponding to two simple 
models for the radiofrequency variation of the absorption factor: 
(a) A = A co- 2 , and (b) A = A (c^-h-p 2 )" 1 . 



or 



N(y)=Aol^8(v-p) 



(4.4) 



This tells one that the consequence of (4.3) is that 
the electron density must be zero everywhere except 
at the height corresponding to the value v=p, 
where N must be infinitely large. This behavoir 
is shown in figure 31), and represents an unrealistic 
solution for the ionosphere. 

C. Example 3. Another model for the absorption 
factor which preserves the asymptotic factor be- 
havior A(u)^(a~ 2 for high frequencies, 



A(o 



(ap 



tan" 



(?) 



(4.5) 



is also plotted in figure 2. The 
this function times co gives: 



Laplace inverse of 



, r/ x 1 dv 1 C° 
v ay* Jo 

_i dv_\ r 

v dy-w Jo 



dx cos (xv) — X 



[--©l 



dx cos (xv)A 



or 



N(v)- 



1 dv_ Aq 
v dy p 



sin xp 
xp 

(l/2,v<p 

l/4,v=p 
0,v>p. 



(4.6) 



The result represents an electron density model 
which is zero below the height corresponding to a 
collision frequency p=v, and constant above this 
height. This N versus v relation is shown in figure 
4 a, and the corresponding N versus y plot for an 
exponential collision (atmosphere) profile v=v 
exp—y/H is given in figure 4b, using the fact that 



v dy" 






On the other hand, this model represents a relatively 
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Figure 4. N versus v relationship corresponding to A(co) = A / 
cop tan -1 (p/co) and equivalent electron density profile for an 
exponentially decreasing collision frequency profile. 



idealized electron density profile and serves to point 
up the profound effect of the low frequency portion 
of the absorption factor A(u>). 

D. Example 4. A fourth example which also 
preserves the a>~ 2 asymptotic behavior 

coyw ~tp 
is plotted in figure 1. The corresponding electron 
density profile is given by 

=- -j — dx cos (xv)A J (xp) 
v ay -k j 



or 



v ay 7r 




1 



Jtf- 






(4.8) 



This model also represents an electron density which 
vanishes below a height corresponding to v=p } but 
has a more interesting variation above there, as 
shown in figure 5a. The corresponding N versus y 
profile for an exponential collision frequency profile 
is shown in figure 5b. The plots again show the 
extreme sensitivity of the electron density profile to 
the detailed behavior of the absorption function at 
the low frequency end. 

5. Alternate Solution 

The fundamental solution (3.8) of the absorption 
integral (3.1) can be simplified in an important way. 
If we again neglect the gyrofrequency terms, we have 

N(y)=- v ~j"dx cos {xv)£ "V1(«)L (5.1) 




Figure 5. N versus v relationship corresponding to A = A / 
coV^ 2 + p 2 and equivalent electron density for an exponential 
collision frequency profile. 

Further reduction is achieved by carrying out the 
cosine transform on the inverse Laplace operation in 
a way which is independent of the function co-A(co). 
To do so, we use the usual complex representation 
for the cosine function and introduce an infinitesimal 
attenuation term in the x integration. 



A7/ , Idv 7 . 1 f : 
vdy <->o 7tJ 

Idv 1 



dxe~ 



v +e- 



ixv\ 



-)x 



cA(w)] x 



lim ( dx[e~ xU ' u) 
v dy 2tt ^o Jo 



+ e 



-xU + iv) 



]£-\vA(co)] x . (5.2) 



Cast in this form, one sees that the integration has 
the effect of undoing the inverse Laplace transform 
and simply to evaluate the function at those values 
given by the coefficients of x in the exponential 
factors, 

N(y)=-ilim^- {M( w )W+M(co)W. (5.3) 
v ay €^o at 

The expression provides an alternative solution to the 
basic absorption (3.1), and could have been derived 
more directly by contour integration techniques. 

The remarkable feature of our alternate solution 
(5.3) is that it depends on the value of the absorption 
function A(oi) only at the complex points e±iv(y), 
where y is the altitude at which Nis to be computed. 
However, one is reminded that the absorption func- 
tion is measured by experiments only for real positive 
values of the carrier frequency w. To exploit the 
new solution (5.3), one must use analytic continua- 
tion from the real positive axis to the complex points 
e±iv, as shown in figure 6. In general, this analytic 
continuation will be determined by all measured 
values of A(u). However, figure 6 suggests that 
the values at e±iv will be especially sensitive to 
values of A(u>) near the origin, and this low frequency 
sensitivity of the inversion process has already been 
noted empirically in the examples of section 4. We 
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Figure 6. Complex frequency plane showing points at which 
absorption function is to be evaluated and emphasizing that 
measured values all lie along positive real axis. 

thus have a more succinct demonstration of the 
limitation of the nondeviative solution considered 
here and suggests the more difficult problem of 
deviative absorption which will recognize the be- 
havior of the absorption-reflection process at low 
frequencies. 

As a concluding exercise, it may be useful to 
resolve the first example of section 4 with the 
alternate solution. If we again take A(o:)=Aqco~ 2 j 
solution (5.3) gives 

vdy e->o 2w\e— iv e + ivj 
1 dv A 



ldvA 0j . ( e \ 
- -j Urn ( 2 , 2 I 



1 dv A «, ~s 



(5.4) 



where we have used here the limit definition of the 
Dirac delta function. Of course, this solution re- 
produces the earlier result (4.2) as one would expect. 



This problem was suggested by some experimental 
work of G. C. Little and his colleagues, which is now 
in progress. Its solution benefited from conversa- 
tions with James R. Wait and M. L. Goldberger. 
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